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Abstract 

Wind tunnels typically have higher 
free stream turbulence levels than are 
found in flight. Turbulence intensity has 
been measured to be 0.5% in the NASA Lewis 
Icing Research Tunnel (IRT) with the cloud 
making sprays off and around 2% with cloud 
making equipment on. Turbulence intensity 
for flight conditions was found to be too 
low to make meaningful measurements 
(<0.1%) for smooth air. This difference 
between free stream and wind tunnel 
conditions has raised questions as to the 
validity of results obtained in the IRT. 
One objective of these tests was to 
determine the effect of free stream 
turbulence on convective heat transfer to a 
smooth and rough airfoil. Another objective 
was to obtain needed heat transfer data for 
the NASA Lewis LEWICE ice growth prediction 
code. 

These tests provide in-flight heat 
transfer data for a NACA-0012 airfoil with 
a 53 3 cm (21 inch) chord. Future tests 
will measure heat transfer from the same 
airfoil in the Lewis Icing Research Tunnel. 
Roughness was obtained by the attachment of 
small, 2 mm diameter hemispheres of uniform 
size to the airfoil in three different 
patterns. Heat transfer measurements were 
recorded in flight on the NASA Lewis Twin 
Otter Icing Research Aircraft. 
Measurements were taken for the smooth and 
roughened surfaces at various aircraft 
speeds and angles of attack up to four 
degrees. Results are presented as Frossling 
number versus position on the airfoil for 
various roughnesses and angles of attack. 


Nomenclature 

A surface area of gage 
c chord length 

d equivalent leading edge diameter 

Fr c Frossling number based on chord 
Fr d Frossling number based on 
equivalent diameter 

h con convective heat transfer coefficient 
k thermal conductivity of air 

M Mach number 

Nu c Nusselt number based on chord 
Nu^ Nusselt number based on equivalent 
diameter 

^Deceased 

*^£raduate Student 

Professor, Chemical Engineering 


Q ei electric power input to heater 
Q en d heat loss from unguarded end of gage 
Q qa n heat loss through gap 
^rad heat loss due to radiation 
Re c Reynolds number based on chord 
Re^ Reynolds number based on equivalent 
diameter 

T t total temperature 
T w measured gage temperature 
v velocity 

o Stephan-Boltzman constant 
/x viscosity of air 
6 density of air 

e surface emissivity of polished 
aluminum 


Introduction 

The hazards of aircraft icing are well 
known. The affect on flow characteristics, 
especially the reduction of maximum lift 
and increase of drag on an iced airfoil are 
documented in references 1-5. A thermal 
analysis of ice accretion shows that 
convective heat transfer is a significant 
factor in the icing process. When 
supercooled water drops from a cloud strike 
an airfoil or engine inlet lip, the heat of 
fusion must be removed before they can turn 
to ice. If the convective, conductive and 
evaporative cooling, as well as the warming 
of impinging droplets, can sufficiently 
overcome the kinetic heating and remove 
enough heat from the water droplets on the 
airfoil surface, then ice will form 6 ' 7 . 
The dominant term in this heat balance is 
convective cooling. 

Therefore icing facilities and ice 
accretion modeling codes must accurately 
reproduce and simulate convective heat 
transfer in natural icing conditions. Wind 
tunnels typically have higher free stream 
turbulence levels than are found in flight. 
Turbulence intensity has been measured to 
be 0.5% in the NASA Lewis Icing Research 
Tunnel (IRT) with cloud making sprays off, 
and around 2% with cloud making equipment 
(nozzle atomizing air only) operating. 
Turbulence intensity for smooth air flight 
conditions measured during this study have 
been found to be too low to make 
meaningful measurements (<0.1%). Somewhat 

higher levels (0. 2-0.4%) were obtained in 
rough air below a layer of cumulus clouds. 
This difference between free stream and 
wind tunnel turbulence has raised questions 
as to the validity of results obtained in 
icing wind tunnels. 


One objective of the present tests is 
to determine the effect of free stream 
turbulence on convective heat transfer to a 
smooth and roughened airfoil. A second 
objective of this work is to obtain much 
needed heat transfer data both for NASA 
Lewis* LEWICE 8 ice growth prediction code 
and to describe the fluid and thermal 
physics occurring during the icing process. 
The NASA code currently uses cylinder in 
crossflow heat transfer data for the 
stagnation region and flat plate heat 
transfer coefficients for the rest of the 
airfoil surface. The present tests provide 
flight test heat transfer data for a NACA- 
0012 airfoil for both a smooth surface and 
three quantifiable roughness patterns. 

Little data presently exists on 
convective heat transfer from an airfoil. 
A NACA study (1946-1951) 9 ' 10 compared in- 
flight convective heat transfer from an 
airfoil, in clear air and during icing 
conditions, with results from the IRT. For 
the flight data two separate airfoils, a 
NACA-0012 and a NACA 65,2-016, were tested 
at zero angle of attack. Only the 65,2-016 
was subsequently tested in the IRT. In 
the "flat plate" region (i.e. the region 
away from the stagnation area) the data 
showed a substantial difference between 
flight and IRT heat transfer on the forward 
portion of the airfoil where the boundary 
layer was laminar. The IRT data was over 
30% higher than the flight data. This 
difference has been attributed to the 
higher turbulence intensities present in 
the IRT. This conclusion is also supported 
by the fact that the flight and IRT data 
agreed fairly well on the downstream 
portion of the airfoil where the boundary 
layer was assumed to be fully turbulent. 

Besides being restricted to a zero 
angle of attack, two other factors limit 
the usefulness of the data for computer 
code predictions. First, the data is 
incomplete and somewhat inconsistent in the 
stagnation region, the area where ice 
growth initiates. Secondly, this data was 
not taken for a rough surface, which can 
significantly alter boundary layer 
characteristics and thus the convective 
heat transfer. Roughness, the result of 
early ice growth, may force a laminar 
boundary layer into transition in the ice 
formation zone. This behavior was observed 
in recent experiments performed on a 
cylinder in crossflow under different 
turbulence and roughness conditions 11 . 
Hence the background turbulence of the IRT 
may not hinder the simulation of airfoil 
ice accretion in flight. 

Heat transfer coefficients on a smooth 
NACA-0012 airfoil in a subsonic wind 
tunnel, as well as on a five minute ice 
accretion shape were measured in reference 
12. The smooth airfoil measurements were 
taken at various angles of attack (-8 
through +8 degrees) and for a chord based 
Reynolds Number range of 7 . 6xl0 5 -2 . 0x10 . 
While the zero degree angle of attack data 
agreed generally with the NACA study, the 
data showed a much larger angle dependence 
on the suction side as compared to the 
pressure side. The data also demonstrated 


a Nusselt Number increase proportional to 
the square root of the Reynolds Number. 

The present study focused on 
convective heat transfer measurements on a 
NACA-0012 airfoil. The NACA-0012 was chosen 
because it is a symmetric profile that is 
commonly used in helicopter main rotor and 
tail rotor applications where ice growth is 
not controlled by electric heating or 
pneumatic boots. Local heat transfer 
coefficients were calculated from 
measurements taken on a smooth and 
roughened NACA-0012 airfoil with a 0.533 
meter (21 inch) chord length. Roughness 
was obtained by the attachment of small 
hemispheres of uniform size (2 mm diameter) 
onto the airfoil in a set and reproducible 
pattern. Three separate position 

patterns, similar to those employed by 
Schlicting 13 in his boundary layer work, 
were used. These patterns were chosen to 
facilitate numerical modeling of the 
roughness in various computer codes. Heat 
transfer measurements were recorded in 
flight on the NASA Lewis Twin Otter Icing 
Research Aircraft. Data were collected 
for smooth and roughened surfaces at 
various aircraft speeds and angles of 
attack up to four degrees. Results are 
presented as Frossling Number versus 
position on the airfoil for various 
roughnesses and angles of attack. Also 
stagnation region data is compared with 
Frossling* s cylinder in crossflow 
solution 14 . For comparison purposes, 

similar tests are also planned in the NASA 
Lewis Icing Research Tunnel. 


Aircraft 

The NACA-0012 airfoil was flown atop 
the NASA-Lewis Twin Otter Icing Research 
Aircraft. The aircraft with the airfoil 
mounted is shown during an aerodynamic 
check flight in figure 1. The Twin Otter is 
a typical twin engine commuter type 
aircraft powered by two 550 shaft 
horsepower turboprop engines. The maximum 
sustainable speed with the NACA-0012 
research airfoil mounted was around 69- 
m/sec (135 knots) at 1585-2250 meters 
(5200-74 00 feet) pressure altitude and a 
temperature range of 289-294 K (60-70°F) . 
The airfoil was mounted on the aircraft by 
attaching the lower end of it to a column 
that extended through the research hatch to 
the floor of the fuselage. The upper end 
of the airfoil was secured by flying wires 
that were attached to the sides of the 
fuselage. 

Airspeed was measured using the 
pitot-static probe built into the boom 
attached to the nose of the aircraft shown 
in figure 1. 

Angles of attack and yaw were measured 
using four pressure sensing ports in the 
hemispherical tip of the boom. The pressure 
difference from the two vertically opposed 
pressure taps was calibrated to measure 
aircraft angle of attack by comparing it to 
deck angle measured with an inclinometer. 
The zero yaw delta-p obtained from the 
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horizontally opposed pressure taps was 
calibrated by aligning a string attached to 
the nose of the aircraft which followed the 
airstream with the aircraft centerline. 
The slope of the yaw delta-p versus yaw 
angle was assumed to be the same as that 
for the angle of attack. 

Free stream static temperature was 
measured with a commercially available 
temperature probe which contained a 
platinum resistance thermometer in a 
specially designed scoop housing. The 
manufacturer supplied calibration data to 
obtain static temperature from the recovery 
temperature measured by the probe and the 
true airspeed. 5 Total temperature was 
calculated using the one-dimensional energy 
equation for a perfect gas under isentropic 
conditions. 

A previous calibration of airspeed 
measured at the boom versus airspeed 
measured at the location of the airfoil was 
used to obtain free stream velocity, total 
temperature and static pressure at the 
test airfoil location. 


Test airfoil 

Heat transfer measurements were made 
on a NACA-0012 airfoil that was designed 
specially for that purpose. The airfoil 
had a chord of 533 cm (21 inches) and a 
span of 1.8 m (6 feet). The airfoil was 
constructed of mahogany and had two spars 
of square, hollow, steel tubing imbedded in 
it. 

An array of heat transfer gages was 
located in a removable section at the 
center of the span. The gages were 
constructed of aluminum and were 6.60 cm 
(2.60 inches) long in the spanwise 
direction, 0.476 cm (0.1875 inch) wide in 
the flow direction, and 0.318 cm (0.125 
inch) deep. Each gage had a groove 
machined into its edge which contained a 
type E (chromel-constantan) , stainless 
steel sheathed, closed, grounded ball 
thermocouple which was held in place with 
an aluminum filled epoxy. A commercially 
available, thin foil heater was fastened to 
the back of each gage with a pressure 
sensitive adhesive. The heat transfer gages 
were held in place with an epoxy that was 
filled with hollow glass microspheres and 
colloidal silica which made final 
contouring to the airfoil profile easier. 
Guard heaters were located beneath the heat 
flux gages to keep heat from leaking out 
the back side of the airfoil. The airfoil 
and epoxy around the gages were sprayed 
with a thin layer of epoxy to seal them 
from moisture. The surface of the gages was 
not coated but was polished to a high 
luster with a polish made for aluminum. 

Figure 2 . shows a cross section of the 
airfoil and the location of the heat 
transfer gages. Table I gives the surface 
distance from the geometric stagnation 
point to the center of each gage and its 
heat transfer surface area. 

The airfoil actually contained 28 heat 
flux gages but only 12 gages in the 


stagnation region were used in these tests 
due to difficulty with the automatic 
controller and data acquisition system. It 
was felt that the gages in the stagnation 
region were of the most interest because 
this is the area where the ice initiates. 
Of these 12 gages only 10 were used to 
report data; gages 1 and 12 were used as 
guard heaters to limit the amount of heat 
leaked from the measuring gages. 

The airfoil was also instrumented with 
two static pressure taps. These taps were 
located on opposite sides of the airfoil at 
the 12% chord position. They were used to 
obtain a measure of angle of attack but 
were not calibrated for that purpose. 

Surface roughness was added to the 
airfoil by fastening hemispheres of silver 
alloy to the surface with cyanoacrylic 
adhesive. The hemispheres were 2 mm in 
diameter and were attached to the airfoil 
in three different patterns. A photograph 
of a typical pattern is shown in figure 3. 
Figure 4 shows sketches of the location of 
the roughness elements relative to the heat 
flux gages for the three patterns. The 
thermal resistance of the gage surface was 
not altered significantly because of the 
sparse spacing of the elements and the high 
conductivity of the silver alloy. No 
attempt was made to account for the 
presence of the roughness elements in the 
data reduction. 


Data acquisition system 

Data collection and recording was 
controlled by a microcomputer. All 
parameters necessary to calculate aircraft 
true airspeed, total temperature, pressure 
altitude, angles of attack and yaw were 
scanned by a commercially available unit 
which contained a multiplexer, signal 
conditioning amplifiers, and a 12-bit 
analog to digital converter. Voltages and 
currents from the heat flux gages were also 
digitized with this unit. Digitized signals 
from this unit were passed to the 
microcomputer which scanned and recorded 
each channel 10 times for each data point. 

Thermocouple extensions were 
terminated at a constant temperature 
reference block whose temperature was read 
with a calibrated platinum resistance 
thermometer. Individual thermocouple 
channels were switched, using a relay type 
multiplexer, to a digital multimeter that 
was capable of reading down to 1 microvolt. 
The IEEE-488 output from this multimeter 
was then recorded by the microcomputer. 
Each thermocouple channel was also scanned 
and recorded 10 times. 


Test procedure 

Turbulence measurements 

Turbulence measurements were made on 
two different flights. The first flight 
was a preliminary test without the test 
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airfoil in place to determine if turbulence 
from any part of the aircraft structure 
would interfere with the heat transfer 
measurements. A constant temperature hot 
wire was mounted about 0.9 meter (3 
feet) above the fuselage in the same 
position as the test airfoil heat flux 
gages. The aircraft was flown during 
daylight both in smooth air and under a 
layer of cumulus clouds. The second 
flight took place in darkness with the 
test airfoil in place. For this test, 
the hot wire was mounted about 2 . 8 
meters (9 feet) forward of the airfoil 
and slightly offset from the aircraft 
centerline. For both flights the hot 
wire was operated in the uncalibrated 
mode as described in reference 16. To 
obtain the bridge voltage at zero 
velocity, the hot wire probe was 
covered with a plastic cylinder and 
allowed to come to equilibrium 
temperature with the air stream. 

At very low turbulence intensities 
(less than about 0.2%), hot wire 
measurements are subject to several 
sources of errors that are not 
important at higher intensities. 
Vibrations of the prongs that the wire is 
mounted on and vibrations of the wire 
itself are among the causes of high 
frequency fluctuating signals that can 
be interpreted as turbulence if one 
only measures the bridge rms voltage 
(spectral analysis of the signal is 
required in order to determine if these 
false signals are present) . To eliminate 
some of the effects of these false signals, 
the bridge voltage was run through a low 
pass filter that was set to cut the signal 
at frequencies above 5 kHz. 

For both flights, the level of 
turbulence intensity in smooth air was 
measured to be around 0.1%. From 
experience with hot wire equipment in 
low turbulence wind tunnels and 
examination of the bridge signal on an 
oscilloscope, it was felt that the 
turbulence intensity for these flights was 
as close to zero as one can get even though 
the numbers from the hot wire equipment 
indicate otherwise. For the flight under 
the layer of cumulus clouds, the intensity 
was measured to be between 0.2 and 0.4%. 
This increased intensity was probably due 
to large scale fluctuations the aircraft 
flew through. It was determined from the 
hot wire measurements that turbulence 
generated by the aircraft structure was not 
a problem and that there was no change in 
the intensity at any of the different 
flight conditions. 


Heat transfer measurements 

All heat transfer data acquisition 
flights were made in darkness to avoid 
solar radiation on the gages and airfoil. 
Flights were conducted at an arbitrary 
altitude that provided smooth atmospheric 
conditions. At low speeds, flaps were 
deployed to minimize the aircraft deck 


angle. At 36 m/sec (70 knots) the measured 
angle of attack was about 1.5 degrees. This 
small angle of attack resulted in a 
slightly swept back test airfoil; this 
effect was ignored in analysis of the data. 

When steady conditions were 
established, the heaters were all adjusted 
to a constant temperature which was 
typically in the range of 306-314 K 
( 90-105 °F) . The heat flux gages were 
operated in the constant temperature 
mode. The temperature of each gage was 
controlled by a circuit that sensed 
thermocouple voltage, amplified it, 
compared that to a reference voltage and 
adjusted the heater voltage to maintain the 
desired temperature. The gain of each 
amplifier could be changed to adjust the 
temperature of individual gages. The 
reference voltage was common to all 
circuits and could be changed to increase 
or decrease the temperature of all gages 
simultaneously. When steady conditions were 
reached, data recording was initiated. 
About two minutes was required to obtain 
and record the required 10 scans of all 
data channels. 

To obtain data for various angles 
of attack on the research airfoil, the 
pilot yawed the aircraft (aircraft yaw = 
research airfoil angle of attack) using a 
combination of rudder and aileron. The 
difference in pressure between the two 
static pressure taps on the airfoil gave a 
measure of the angle of attack; this 
quantity was recorded with the other data. 
Figure 5 is a plot of aircraft yaw angle 
measured from the boom delta-p versus the 
pressure difference between the two static 
taps on the airfoil made dimensionless by 
dividing by the dynamic pressure (<$v 2 /2) . 
This plot gives a measure of the scatter in 
the angle of attack. The scatter in the 
abscissa is the accuracy with which the 
pilot could set and hold the aircraft yaw. 


Data reduction 


The average heat transfer coefficient 
from each gage was obtained from the 
applied heater voltage and current, and the 
temperature difference between the gage 
and the calculated free stream total 
temperature. Since only the convective 
heat transfer was desired, the radiation 
heat loss had to be subtracted from the 
total electric power input to each heater. 
Further, the heater gages embedded in the 
airfoil were secured in place and separated 
from each other by an epoxy resin. Some 
heat was conducted from the edges of each 
gage, through the epoxy and convected from 
the surface of the foil in the gaps between 
gages and from the unguarded ends of the 

gages. These losses were also subtracted 
from the electric power. Therefore, the 
local convective heat transfer coefficient 

for each aluminum heater gage was 

calculated from: 

Qe I -Q r ad~Qgap“Qend 

n con “ ( x ) 


A(T W - T t ) 
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where Q EI (voltage X current) is the total 
electric power input to each heater. The 
quantity Q rac j is the radiation heat loss, 
which is estimated by: 

Qrad = oAe(T w 4 - T t 4 ) (2) 


A value of 0.045 was used for e, the 
emissivity of polished aluminum. Q c 

and Q end are the heat loses through 
epoxy gaps separating the aluminum gages 
and the unguarded ends of the heaters. 
These were obtained from an exact solution 
for heat conduction in a rectangle with 
appropriate boundary conditions as 

detailed in reference 11. The remaining 
quantities are: A, the surface area of 

each aluminum gage, T„, the measured gage 
temperature, and T t , the total temperature 
calculated from the measured static 
temperature and the true airspeed, i.e. 


T t = T s (1+M 2 /5) (3) 


where M is the Mach Number. Two Frossling 
Numbers were employed in this analysis, one 
based on chord length and the other based 
on an equivalent leading edge diameter. 
This equivalent diameter is defined as the 
diameter of a cylinder inscribed in the 
leading edge of the airfoil. The Frossling 
Number based on chord was calculated as: 

Nu c ( h con c / k > 

Fr = = (4) 

7Re c 7(5v c/n) 

where c is the 0.533 meter (21 inch) chord 
length. The Frossling Number based on 
equivalent diameter was calculated as: 


Nu d (^con^/^1 , _ 

Fr d = = (5) 

7Re d 7(6vd//i) 


where the equivalent diameter, d, for a 
NACA-0012 airfoil is 3.16% of the chord 17 
or 1.69 cm (0.664 inch) for the airfoil 
tested. The density, 6 , was calculated 
from the ideal gas relation for air using 
the static temperature and pressure at the 
test airfoil location. Velocity was the 
calculated true airspeed at the test 
airfoil location. The thermal conductivity, 
k, and viscosity, fi, were obtained as 
functions of temperature from the air data 
of reference 18. These thermal properties 
were evaluated at the film temperature 
given by: 

T f = (T w + T t )/2. (6) 


Error analysis 

An error analysis according^ to the 
method of Kline and McClintock 19 was 
performed on each of the calculated local 
convective heat transfer coefficients, 
hcon • The e ^ rors for each gage were 
similar and averaged around 15%. The 
majority of this error was found to be due 
to uncertainty in the gap heat loss term 
because the thermal conductivity of the 


epoxy gaps was not known and had to be 
estimated. This would not be a random 
error but would tend to bias the data 
either high or low. 


Results and Discussion 

In this section heat transfer data for 
smooth and roughened airfoil surfaces will 
be presented as Frossling number based on 
chord length versus dimensionless surface 
distance from stagnation point (s/c) . 
These results will be presented for nominal 
0, 2, and 4 degree angles of attack, and 
for various Reynolds numbers in the range 
1.24 x 10 6 to 2.50 x 10 6 . Table 2 
contains the Frossling numbers for all 
gages, Reynolds numbers, and angle of 
attack for data points taken. A comparison 
will be made of the present data with 
previous airfoil heat transfer work and 
with an analytical solution for a cylinder 
in cross flow. 

Smooth airfoil 

Figure 6 shows Frossling number based 
on chord as a function of s/c for the 
smooth airfoil at 0 degree angle of attack 
for several Reynolds numbers. The data 
plotted in this manner collapses onto a 
single curve which shows that the heat 
transfer coefficients are proportional to 
the square root of Reynolds number. The 
solid line on the figure represents the 
averaged, smooth-surface, zero degree angle 
of attack data and will be plotted on 
subsequent figures for reference. As 
expected the Frossling number is greatest 
at the stagnation point, an average value 
of 4.3 being observed, and then trailing 
off smoothly to an average value of 1.0 at 
s/c of 8.3. The "bump* 1 at s/c = 4.8 cannot 
be explained; there are no obvious 
roughnesses or steps in the surface at this 
or any other point. Perhaps there is a 
subtle anomaly in the profile at this point 
that has not been detected. 

Defining Frossling number in terms of 
an airfoil leading edge equivalent diameter 
allows comparison of Frossling* s analytical 
solution for heat transfer in the 
stagnation region of a circular cylinder 14 
with the present data. The average 
Frossling number based on leading edge 
equivalent diameter for the smooth airfoil 
was found to be 0.76, roughly 22% lower 
than the 0.945 value predicted by 
Frossling* s cylinder solution. Frossling *s 
analytical results are often used with an 
equivalent leading edge diameter to compute 
heat transfer in the stagnation region for 
airfoils and turbine blades but, to our 
knowledge, no one has ever conducted an 
experiment to prove the validity of this 
method . 

Figure 7 is a plot of Frossling number 
based on chord against s/c for the smooth 
airfoil at a 2 degree angle of attack. Data 
for the suction side of the airfoil are 
represented by the positive s/c values. 
This convention will be maintained 
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throughout this report. The data again 
collapses onto a single curve and 
illustrates good agreement with the yRe c 
dependence. Comparison with the 0 degree 
smooth foil data shows no notable 
difference. 

Figure 8 shows Frossling number for a 
smooth airfoil at 4 degrees angle of 
attack. Again the data for all Reynolds 
numbers can be represented by a single 
curve. Comparing with 0 and 2 degree data, 
shows very little angle dependence save a 
slight increase (11%) on the leading edge 
and a slight decrease (6%) on the first 
gage of the pressure side of the airfoil. 
This behavior can be explained by the 
movement of the aerodynamic stagnation 
point toward the pressure side with 
increasing angle of attack. This changes 
the radius of curvature at the aerodynamic 
stagnation point; the stagnation region 
sees an effectively larger equivalent 
diameter and this results in a lower heat 
transfer coefficient. The flow is then 
highly accelerated around the leading edge 
increasing heat transfer at the geometric 
stagnation point. 

Leading edge roughness 

Figures 9,10, and 11 show Frossling 
number based on chord versus s/c for an 
airfoil with roughness elements attached to 
the leading edge, as shown in figure 4b.), 
for angles of attack equal to 0, 2, and 4 
degrees, respectively. The roughness 
element row positions are denoted by the 
arrows below the abscissa. Compared with 
data from the smooth airfoil, the leading 
edge roughness increases the heat transfer 
an average of 8% at the stagnation point 
but it returns to essentially the smooth- 
surface values away from the stagnation 
region. This increase may be partially 
explained by the 4% increase in surface 
area due to the presence of the 
hemispherical roughness elements. It could 
also be attributed to a disturbance of the 
boundary layer by the relatively huge 
roughness elements followed by a return to 
laminar flow sufficiently past the leading 
edge trip point. Note also that this data 
set exhibits the same slight angle of 
attack dependence for the first pressure 
side and leading edge gages as the smooth 
airfoil. The first gage on the suction side 
measures slightly (8%) higher heat transfer 
for the 4 degree angle of attack than for 0 
or 2 degrees. 

Sparse roughness 

Frossling number as a function of s/c 
for the sparse roughness pattern shown in 
figure 4 c.) at 0, 2 and 4 degrees angle of 
attack are presented in figures 12, 13, and 
14 respectively. Compared with the 
smooth airfoil, this data set contains 
somewhat more scatter but still can be 
represented by a single curve which is 
consistent with the 7Re p dependence. For 
the 0 degree angle of attack case, the heat 
transfer increase at stagnation (gage 4) is 
11%; comparable the leading edge roughness 


case. Past the stagnation region, the heat 
transfer exhibits a pattern of increasing 
at and immediately downstream of the 
roughness row position and then falling off 
slightly. At gage 7, the heat transfer 
increases by 55% over the smooth airfoil 
case, increases to 59% at gage 8, then 
falls to a 52% increase at gage 9. At gage 
10 the next row of roughness elements is 
encountered. The heat transfer at gages 10 
and 11 increases by about 170% over the no 
roughness case. The sensitivity of the 
boundary layer to roughness seems to 
increase with downstream location. 

The angle of attack dependence is also 
more pronounced for the sparse roughness 
pattern than the smooth and leading edge 
roughened cases. For the sparsely 
roughened airfoil, the Frossling numbers 
for 2 degrees angle of attack for gages 7 
through 11 increase gradually with s/c from 
8 to 15% over the 0 degree case and from 14 
to 26% going from 0 to 4 degrees. 
Increasing angle of attack causes heat 
transfer to increase with s/c over the 0 
degree case. Note that for 4 degrees, the 
characteristic increase at stagnation and 
slight decrease on the second gage of the 
pressure side of the airfoil are also 
observed. 

Dense roughness 

Frossling number versus s/c for the 
dense roughness pattern (fig. 4 d.)) at 0, 
2 and 4 degrees angle of attack are 
presented in figures 15, 16, and 17. For 
the 0 degree angle of attack case, the data 
points still tend to fall on one curve 
indicating a 7Re c dependence. Increasing 
the density of roughness elements from the 
sparse to dense pattern had a dramatic 
effect on heat transfer downstream. For 
the 0 degree angle of attack, gage 6 
increased 3 2 % and gages 7 and 8 increased 
an average of 54% over the sparse roughness 
case. Further downstream past gage 7, the 
density of roughness elements decreases and 
at gages 10 and 11, the effect of the 
increased density of the roughness elements 
upstream seems to have nearly damped out. 
This trend indicates that if there is 
roughness of sufficient magnitude present, 
the boundary layer is perturbed locally and 
immediately downstream but, as the density 
of roughness is lowered in the downstream 
direction, the heat transfer recovers to a 
level that is consistent with the sparse 
roughness pattern. 

For the 2 and 4 degree angles of 
attack there is considerably more scatter 
in the data than was present in the smooth 
airfoil cases. At high values of s/c, the 
Frossling numbers increase monotonically 
with Reynolds number; this may indicate a 
trend away from the ,/Re c dependence with 
increasing roughness and angle of attack. 

The angle of attack dependence is much 
more prominent in the dense roughness case 
compared to other cases tested. A increase 
from 0 to 2 degrees caused roughly a 2 0% 
increase in Frossling number for the gages 
between s/c of 2 and 5, while a 4 degree 
change yielded an increase of roughly 
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39%. For gages at s/c locations greater 
than 5, increasing angle of attack from 0 
to 2 degrees caused a 15.3% increase and 
from 0 to 4 degrees a 27% increase in 
Frossling number. 

Comparison with other data 

Finally figures 18 and 19 compare the 
present smooth airfoil results with 
previously published data. Figure 18 shows 
the comparison with flight test data for 
both a NACA-0012 and a NACA 65,2-016 
airfoils 8 9 . Aside from a few exceptional 
points good general agreement is observed, 
most values agreeing within 10 %. 

Figure 19 illustrates the comparison 
of the present smooth airfoil, 0 degree 
angle of attack flight data with wind 
tunnel data of reference 12. Relatively 
good agreement exists up to a s/c value of 
about 5? however, further down the airfoil 
Frossling numbers differ by nearly 300%. 
This could be due to wind tunnel turbulence 
or roughness of the reference 12 model 
surface. In a personal communication with 
the author of reference 12 it was indicated 
that the surface of their airfoil was 
"rough" . 
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Table 1 . Location and surface area of heat transfer gages 


GAGE# 

S/C 

SURFACE 

AREA 

sq.cm. 

1 

— 0.036 

3.145 

2 

-0.024 

3.145 

3 

-0.012 

3.145 

4 

0 

3.187 

5 

0.012 

3.145 

6 

0.024 

3.145 

7 

0.036 

3.145 

8 

0.048 

3.145 

9 

0.060 

3.145 

10 

0.072 

3.145 

11 

0.083 

3.145 

12 

0.095 

3.145 


Table 2. Frossling number for each gage for all conditions tested. 


roughness 

pattern 

alpha 

Re 

2 

3 

4 

Gage Number 
5 6 

7 

8 

9 

10 

11 

none 

.027 

1291850 

3.140 

3.959 

4.495 

3.960 

2.899 

2.312 

2.003 

1.167 

1.034 

0.988 


-.441 

1924300 

3.027 

3.803 

4.261 

3.787 

2.780 

2.207 

1.915 

1.133 

1.004 

0.954 

none 

-.435 

1935420 

3.077 

3.858 

4.332 

3.867 

2.806 

2.208 

1.929 

1.145 

1.003 

0.954 

none 

-.127 

1935620 

2.935 

3.691 

4.174 

3.703 

2.716 

2.169 

1.866 

1.122 

0.971 

0.919 

none 

-.310 

2482730 

3.016 

3.752 

4.190 

3.798 

2.804 

2.225 

1.938 

1.154 

0.999 

0.925 

none 

-2.60 

1285880 

3.151 

3.949 

4.646 

3.906 

2.825 

2.239 

1.892 

1.151 

0.993 

0.958 

none 

-2.33 

1937560 

3.079 

3.845 

4.410 

3.907 

2.694 

2.140 

1.861 

1.111 

0.956 

0.886 

none 

-2.46 

2482670 

3.061 

3.794 

4,289 

3.761 

2.715 

2.150 

1.816 

1.046 

0.930 

0.873 

none 

-4.41 

1273580 

3.210 

3.747 

4.916 

3.943 

2.765 

2.157 

1.826 

1.097 

0.948 

0.899 

none 

-4.49 

1283910 

3.189 

3.734 

4.898 

3.912 

2.721 

2.107 

1.784 

1.075 

0.915 

0.866 

none 

-4.35 

1934040 

3.106 

3.504 

4.749 

3.727 

2.620 

2.071 

1.736 

1.031 

0.903 

0.824 

none 

-4.32 

2465000 

2.966 

3.315 

4.505 

3.719 

2.622 

2.064 

1.733 

1.056 

0.898 

0.846 

l.e. 

-.568 

1258000 

3.165 

4.069 

4.857 

4.065 

2.938 

2.341 

2.002 

1.181 

1.045 

0.996 

l.e. 

-.266 

1886310 

3.104 

3.921 

4.662 

3.866 

2.817 

2.251 

1.921 

1.168 

1.004 

0.927 

l.e. 

-.423 

2424730 

3.156 

3.973 

4.730 

4.014 

2.907 

2.313 

1.898 

1.167 

1.015 

0.949 

l.e. 

-2.32 

1265320 

3.177 

4.018 

4.872 

3.973 

2.850 

2.246 

1.933 

1.153 

1.009 

0.952 


-2.34 

1878440 

3.005 

3.771 

4.577 

3.725 

2.692 

2.194 

1.848 

1.078 

0.966 

0.866 

l.e. 

-2.42 

1884740 

3.031 

3.833 

4.571 

3.766 

2.700 

2.124 

1.830 

1.072 

0.963 

0.878 

l.e. 

-2.32 

2392890 

3.028 

3.825 

4.519 

3.881 

2.771 

2.443 

1.814 

1.090 

0.975 

1.142 

l.e. 

-2.31 

2396510 

3.039 

3.817 

4.520 

3.850 

2.867 

2.321 

1.864 

1.039 

0.946 

0.836 

l.e. 

-4.78 

1244970 

3.155 

3.649 

5.158 

4.064 

2.819 

2.181 

1.863 

1.128 

1.008 

1.003 

l.e. 

-4.33 

1879440 

2.998 

3.442 

4.933 

4.208 

2.655 

2.016 

1.478 

1.207 

0.953 

0.898 

l.e. 

-4.37 

24224 80 

3.066 

3.506 

5.089 

4.131 

2.860 

2.203 

1.859 

1.166 

1.115 

1.216 

sparse 

-.272 

1257300 

3.254 

4.160 

5.044 

4.205 

2.989 

3.487 

3.071 

1.728 

2.577 

2.502 

sparse 

-.138 

1884970 

3.084 

3.973 

4.736 

3.904 

2.817 

3.406 

3.057 

1.727 

2.690 

2.663 

sparse 

-.333 

1886120 

3.028 

3.911 

4.626 

3.904 

2.714 

3.384 

3.016 

1.715 

2.658 

2.626 

sparse 

-.304 

2415640 

3.085 

3.920 

4.662 

3.936 

2.807 

3.487 

3.142 

1.774 

2.837 

2.816 

sparse 

-1.94 

1260160 

3.299 

4.181 

5.052 

4.129 

2.882 

3.602 

3.173 

1.819 

2.796 

2.765 

sparse 

-2.32 

1896200 

3.147 

3.897 

4.701 

3.911 

2.745 

3.631 

3.276 

1.871 

2.978 

2.946 

sparse 

-2.28 

2420490 

3.110 

3.905 

4.685 

3.970 

2.798 

3.771 

3.444 

1.991 

3.237 

3.232 

sparse 

-2.33 

2461170 

3.156 

3.956 

4.738 

4.041 

2.883 

3.901 

3.565 

2.057 

3.315 

3.325 

sparse 

-4.24 

1244580 

3.417 

3.971 

5.586 

4.433 

3.001 

4.002 

3.591 

2.142 

3.342 

3.177 

sparse 

-4.46 

1263380 

3.342 

3.914 

5.499 

4.322 

2.935 

3.865 

3.473 

2.049 

3.218 

3.062 

sparse 

-4.44 

1890090 

3.122 

3.594 

5.133 

4.099 

2.749 

3.866 

3.579 

2.140 

3.503 

3.294 

sparse 

-4.38 

2434720 

3.001 

3.431 

4.987 

4.103 

2.745 

4.000 

3.783 

2.333 

3.526 

3.466 

dense 

.100 

1269450 

3.732 

4.339 

5.072 

4.737 

3.562 

5.424 

4.668 

2.456 

2.590 

2.439 

dense 

.885 

1276420 

3.771 

4.328 

5.066 

4.708 

3.576 

5.497 

4.679 

2.484 

2.614 

2.427 

dense 

-.318 

1887710 

4.010 

4.081 

4.717 

4.421 

3.814 

5.605 

4.740 

2.498 

2.738 

2.572 

dense 

-.630 

1908240 

3.963 

3.985 

4.696 

4.361 

3.819 

5.564 

4.717 

2.514 

2.710 

2.532 

dense 

-.335 

2446790 

4.136 

3.962 

4.434 

4.261 

3.891 

5.555 

4.761 

2.562 

2.811 

2.657 

dense 

-2.80 

1275150 

3.430 

4.158 

5.026 

4.840 

4.354 

6.295 

5.204 

2.764 

2.891 

2.684 

dense 

-2.21 

1929750 

3.543 

4.011 

4.738 

4.759 

4.549 

6.535 

5.402 

2.896 

3.099 

2.899 

dense 

-2.34 

2466420 

3.660 

3.952 

4.661 

4.809 

4.821 

6.748 

5.721 

3.124 

3.276 

3.064 

dense 

-4.83 

1275190 

3.308 

3.740 

5.379 

5.414 

5.138 

7.273 

5.730 

3.030 

3.157 

2.932 

dense 

-4.24 

1899520 

3.209 

3.566 

5.215 

5.452 

5.473 

7.356 

5.958 

3.230 

3.383 

3.125 

dense 

-4.40 

2439820 

3.283 

3.545 

5.284 

5.699 

5.889 

7.924 

6.460 

3.539 

3.657 

3.381 


8 






I III 


Figure! NASA Lewis Icing Research Aircraft with NACA- 001 2 heat transfer research airfoil 
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Figure 2. Cross section of NACA - 0012 airfoil with heat flux gages 


Figure3. Sparse roughness pattern 
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d.) Dense roughness pattern. 
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b.) Leading edge roughness pattern. 



a.) Typical roughness element. 

Figure 4. Location of roughness elements relative to 
heat flux gages. 



Figure 5. Comparison on angle of attack 
measurements from airfoil and 
aircraft instruments. 
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DIMENSIONLESS SURFACE DISTANCE 

Figure 6. Frossling number based on chord for smooth airfoil, 

zero degree angle of attack for various Reynolds numbers. 



Figure 7. Frossling number based on chord for smooth airfoil, 

two degree angle of attack for various Reynolds numbers. 



Figure 8. Frossling number based on chord for smooth airfoil, 

four degree angle of attack for various Reynolds numbers. 
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I 

Figure 9. Frossling number based on chord for leading edge roughened 
airfoil, zero degree angle of attack for various Reynolds numbers. 



I 

Figure 1 0. Frossling number based on chord for leading edge rougnened 
airfoil, two degree angle of attack for various Reynolds numbers. 



t 

Figure 1 1 . Frossling number based on chord for leading edge rougnened 
airfoil, four degree angle of attack for various Reynolds numbers. 
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Figure 1 2. Frossling number based on chord for sparsely roughened airfoil, 
zero degree angle of attack for various Reynolds numbers. 



I 1 I 

Figure 13. Frossling number based on chord for sparsely roughened airfoil, 
two degree angle of attack for various Reynolds numbers. 



I I 1 

Figure 1 4. Frossling number based on chord for sparsely reoughened airfoil, 
four degree angle of attack for various Reynolds numbers. 
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fit! L 

Figure 1 5. Frossling number based on chord for densely roughened airfoil, 
zero degree angle of attack for various Reynolds numbers. 
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Figure 16. Frossling number based on chord for densely roughened airfoil, 
two degree angle of attack for various Reynolds numbers. 



I I t 1 I 

Figure 17. Frossling number based on chord for densely roughened airfoil, 
four degree angle of attack for various Reynolds numbers. 
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DIMENSIONLESS SURFACE DISTANCE, S/C 


Figure 18. Comparison of NASA 1987 smooth airfoil data with flight test 
data of reference 9. 



DIMENSIONLESS SURFACE DISTANCE, S/C 


Figure 19. Comparison of NASA 1987 smooth airfoil data with wind tunnel 
data of reference 12. 
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